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Abstract 



For a one-parameter family of lower triangular matrices with entries involving Gegenbauer 
polynomials an explicit inverse is given, again with entries involving Gegenbauer polynomials. 
One choice of the parameter solves an open problem in a recent paper by Koelink, van 
Pruijssen & P. Roman. Another family of pairs of mutually inverse lower triangular matrices 
with entries involving Gegenbauer polynomials, unrelated to the family just mentioned, is 
implied by a paper by J. W. Brown & S. M. Roman (1981). J. Koekoek and R. Koekoek 
(1999) generalize this family to entries involving Jacobi polynomials. The present paper 
also shows that this last family is a limit case of a pair of connection relations between 
Askey- Wilson parameters having one of their four parameter in common. 

1 Introduction 

The first part of this note is intended as a kind of supplement to the paper [9j by Koelink, 
van Pruijssen & Roman. It solves their open problem [9l Theorem 2.1 and paragraph after 
Theorem 6.2] to invert a lower triangular matrix with entries involving Gegenbauer polynomials. 
For a one-parameter family of such matrices I give the explicit inverse matrix in Theorem I2.1[ 
One specialization of the parameter gives the inversion desired in [9]. Another specialization 
gives a matrix inversion already handled in a paper [2] by Brega & Cagliero. 

Another two-parameter family of pairs of mutually inverse lower triangular matrices with 
entries involving Gegenbauer polynomials, unrelated to the family just mentioned, is implied 
by J. W. Brown & S. M. Roman O (4.14)]. J. Koekoek and R. Koekoek 13 (17)], unaware 
of [3], generalized a one-parameter subfamily of this two-parameter family to entries involving 
Jacobi polynomials. I will show that this last family can be realized as a limit case of a pair 
of connection relations between Askey- Wilson parameters having one of their four parameter in 
common. These Askey- Wilson connection coefficients were first given by Askey & Wilson [H 
(6.5)]. The limit case connects Jacobi polynomials P^'^^ with shifted monomials x i— >■ (x — y)^ . 

The contents of the paper are as follows. The main results are stated in Section [2j The 
computations leading to the explicit inverse matrix of the first family of lower triangular matrices 
are given in Section [3l The computations giving the limit of the Askey- Wilson connection 
relations are done in Section 21 

Acknowledgement I thank Michael Schlosser for the suggestion to look for a limit case of the 
Askey- Wilson connection relations, and I thank Roelof Koekoek for calling my attention to [8]. 
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2 Main results 



Jacobi polynomials (see for instance |8l Section 9.8]) can be expressed in terms of the Gauss 
hypergeometric function by 

pt'H-) ■■= .i^i f " + ! t ^ ^ - (2.1) 

E (n + a + /3 + l)fc {a + k)n-k jx - 1)'' , . 

k\(n-k)\ 2^ ' ^ ' ' 

k=o ^ ' 

Note that they are well-defined for all values of a, /?. Their normalization avoids artificial 
singularities. For a = /3 Jacobi polynomials are often written as Gegenbauer polynomials: 

C7W(,):= J^Pri^-i)(,). (2.3) 

Thus (x) = 6nfi, which will be kept as a convention in this paper, although in literature the 
case A = is usually rescaled in order to obtain the Chebyshev polynomials of the first kind. 
In Section [3] it will be shown that 

n 

p(a,+kA+k)^^^ p(a,-k,P,^k) ^ p^ai+a2,/3i+/32)(^)_ (2.4) 

fc=0 

and 

( t. n(a2 + P2) \ (a^+k,h+k), . p(a2-fe,/32-fc). X 



fc=0 



ai +a2 + /3i +/32 + 2n ""^ ^ ' 



For ai = /?!, a2 = 1^2 formula ()2.5p reduces to 



±{k- ^^f^) Pi-+'='-+^')(x) p(--'='--'^)(x) = (n > 0). (2.6) 



fe=0 

A further reduction to the case ai = a, 02 = — o (a 7^ 0) yields 

^ + (Q+A;,a+fc)/ N jy(-a-k-a-k) I 



^ Ptk \^) '""^^(^) = ^n,o. (2.7) 



k=0 



Now make in (12. 7p the substitutions n^m — n, k^k — n, a^a + n. The resulting 
identity for a > is 

E ^ n^:r''-"-'^(-) = W (- > n > 0). (2.8) 

fe=n 

This identity can be rephrased as a pair of inverse relations or equivalently as two lower triangular 
00 X 00 matrices with explicit entries which are inverse to each other: 
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Theorem 2.1. AB = I = BA where A,B are lower triangular matrices given by (a > 0) 

Amu = Pln^r^''\xl = ^ Pi-r'"'— "^^(x) {m>n> 0). (2.9) 

n + a 



Just as we can go from AB = I to (j2.8p and backwards, we can go back and forth from 
BA = I to the identity 

n 

S ^ ^.^""^(-) ^i:r'-"-"^(-) = Ko. (2.10) 

k=0 

There are two different places in hterature where Theorem 12.11 can be used, for a = 1 and 
Q = |, respectively: 

1. The case a = 1 occurs in Brega & Cagliero [2, p. 471] with a proof similar as given here. 

2. The case a = ^ of the matrix A in (j2.9p occurs in Koelink, van Pruijssen & Roman [9l 
Theorem 2.1] in the form of the lower triangular matrix L{x) given by 

There the matrix has finite size (which does not matter for the purpose of inversion). As 
the authors wrote in [9l paragraph after Theorem 6.2], they tried to find an explicit inverse 
matrix but did not succeed. We can give the inverse by ()2.9p for a = ^ as follows. 

(2.12) 

Remark 2.2. Brown & Roman [3^, (4.12)-(4.15)] obtain inverse relations involving Gegenbauer 
polynomials of which a special case is close to (j2.7p but not equal to it. It reads 

^ (2a + 2),, W-dn,o- l^.i^j 

In fact, they give a more general identity 



n 



E ^r^^^(-) Ctt-'\x) = 0. (2.14) 

A;=0 ^ 



Then f^HM is the case fi = -1, v = -a - ^ of (ITT^ . 
Formula (|2.13p is also the case q = /3 of the identity 



n 



=*»,0. (2.15) 

This last identity is a consequence of the pair of mutually inverse lower triangular matrices 
implied by J. Koekoek & R. Koekoek [7, (17)]. 

In sectional will show that (j2.15p . and hence ()2.13p . is related to a limit case of a connection 
formula for Askey- Wilson polynomials. 
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Remark 2.3. There remain several interesting questions. First of all, is there a larger family of 
explicit mutually inverse lower triangular matrices which includes both the family of Theorem 
12. H and the family (j4.1Up implying (|2.15|) ? Furthermore, are there two simple systems of special 
functions connected by the matrices in Theorem 12. II .'' If yes, can this also be seen as a limit case 
for g ^ 1 of some connection formula in the g-case? Finally there is the puzzling Brown-Roman 
formula (|2.15p . Does this have an extension to Jacobi polynomials for general /i? 

3 Computations leading to Theorem 12.11 

Lemma 3.1. If the functions f and g have derivatives up to order n then 

E = ifaY'^Hx), (3.1) 

fc=0 ^ ^ 

= n(/5')^"-'H^)- (3.2) 

fc=o ^ ^ 

Proof Formula (j3.ip is well-known. For the proof of (j3.2p rewrite its left-hand side as 

n-l 

n 

j=0 

and use (O)- □ 



E("T')/("-^-^)(x)(<7')(^-)(x) 



Jacobi polynomials ()2.ip satisfy a Rodrigues' formula 

pt^\x) = ^-^{i-xr^{i+xr^ (^)" ((i-^)"+"(i+^)"'''')- (3-3) 



Hence 



^pK+Mi+fc)(^)p^a2-fc,/32-fc)(^) = ^^(l-x)-°i-"2 (l+x)-^i-^2 
A;=0 ^' 



k=0 



By (KTlj) and again (jO]) we obtain (12:4]) 

Similarly, by (j3.3p and (|3.2p we can write for n > 0: 



k=0 ^' 

2*^(71 — 1)' \dx J V ^ 
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By straightforward computation we get 



(1 - (1 + x)"+^i ^ ((^ " (1 + ^)^' 



— (1 -x)"i+"2+n(^^^)^ 



ai + 02 + /3i + /32 + 2n 

ai + 02 + /3i + /32 + 2n ^ 

By (13. 3p we finally obtain 

+ "2/31 - ai/?2 + n(a2 - /?2) ^(ai+aaA+fc) . n . ^ qn 

«! + a2 + /3i + /32 + 2n """^ 



Combination of ([23]) and ([331) yields (I23]) . 



4 Limits of a connection formula for Askey- Wilson polynomials 

Askey- Wilson polynomials [1] are defined by 

(0102,0103, 0104; g)„ /Q'"",ala203a4Q'""^ale*^,ale"*^ \ 

p„(cos 6^; ai, 02, 03, 04 | g) := 4<p3 ; Q • 

aj^^ \ 0102,0103,0104 / 

(4.1) 

They are symmetric in 01,02,03,04. The connection coefficients Cn,k in 

n 

Pn{cos0;bi,b2,b3,a4, \ q) = ^ c„,fc(6i, 62, 63, 04; oi, 02, 03, 04 | g) pfc(cos 61; oi, 02, 03, 04 | q) (4.2) 

k=0 

are explicitly given in Askey & Wilson [H (6.5)]: 

qk{k-n)(^q.q^^ (616263049""^; g)fc 



Cn,k{bi, &2, ^3, 04; oi, 02, 03, 04 I q) 



al ^{q;q)n~k{q]q)k (01020304?^^ 



n k r k u k ^ : (1 ",61 62^3049^ l,Oi04g'',0204g'=,0304g'' \ 

X (6104^ ,62049 ,63049 ;g)n-fc 5(/'4 , , . , . nr. • (4.3) 

\ 61049'', 62049'^, 63049'^, 010203049^'^ y 

See also Ismail & Zhang [5l Section 3] and Ismail [6l §16.4], where the connection coefficients are 
given more generally for 04 7^ 64. However, note that in [5, (3.13)] and [6, (16.4.3)] one should 
read c„^fc(b, a) instead of c„^fe(a, b). 
Now put 

04:=9"+V&i, 63:= 9^+7^2 (4.4) 
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in (112]) and (03]), and multiply both sides of (gj]) by l/{q;q)n. By (gJl and we see that 

lim-^p„(cos0;6i,62,/+V&2,g"+V&i I q) 

(l-6i62)(62-foi)\" (a + l)„ ^ /-n,n + a + /3 + 1 62(1 - 26i cos + 

2-^*1 



6162 / n! V a + 1 ' (1 - &i62)(fe2 - ^i) 

^ / (l-blfe2)(fe2-fcl) \ " («,^) / 62(1 - 25i COS g + 

V ^1^*2 y " V {i-bib,){b2-b,) )■ 

By ()4.2p we also see that 

limpfc(cos6';ai,a2,a3,g°+V&i I 

(1 - aia2)(l - aia3)(6i - ai)\ ( -k {hi - aia2az){l - 2ai cos6' + af) 



aihi ) ^ - ' (1 _ a^a2)(l - aia3)(6i - ai) 

(1 — a\a2)(\ — a\a-i)(})\ — a\) — (61 — aia2a3)(l — 2a\ cos + af ) ^ ^ 



For the 5(^)4 in (j4.3|) we get 

( a+fe+i N _^ /g'=-",g"+'=+°+'3+i,g°+^'+iai/6i,g"+^+ia2/6i,g°+'=+^a3/6i 
hm ; ^ 5(/)4 f 



g^i (g;g)n-fc V <7"+'^+i,g"+^+i62/fei,<7"+'^+'=+V(^i^2),g"+2fc+ia^a2a3/5^ 

(a + /c + l)„„fc ^/-n + A;,n + A; + a + /3 + l^ &2(&i - - a2)(&i - 03) \ 

(n - A:)! a + fc + l ' (61 - 62)(6i&2 - l)(6i - aia2«3) / 

p(a+fc,/3+fc) A „ ^2(^1 - ai)(&i - a2)(&i - 03) 



(61 - h2)(h\h2 - l)(6i - 010203) / ' 
For the other factors in (|4.3p we get 



</^ior''((?;(?)fc ('?°+'=oia203/6i;g)/^ ''^ ' 

(n + a + /3 + l)fc / 61 \V(^i-^2)(&i62-l)^"-*^ 



/c! V^l ~ <^l<^203/ V ^1^2 

Also put 

^ ^ ^ _ 2 62(1-261 cosg + 6f) y ^ 1 _ 2 62(61 -oi)(6i -02)(6i -03) 

(1 - 6i62)(62 - 61) ' (61 - 62)(6i62 - l)(6i - 010203) ■ 

Then we obtain the following limit case of (j4.'2p as 1: 

ft-<"(^) = E (^) ' . (4.6) 

fc=o ■ ^ ^ 
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Now interchange the a and b parameters in (j4.2p : 

n 

Pn{cos9;ai,a2,a3,ai \ q) = ^ c„,fe(oi, 02, as, 04; 61, 62, 63, 04 | g) pfc(cos 61; 61, 62, 63, 04 | q), 

(4.7) 

and use ()4.3p with the a and 6 parameters interchanged and with the order of summation 
reversion formula [T, Exercise 1.4(ii)] appUed to the ^(j)^: 

Cn,fc(ai,a2,a3,a4;6i,62,03,a4 | g) = ^^3^- 



(61049 ,62049 ,03049 ;9)n-fc 



(61626304^^ 1; q)k{bib2bsa4,q'^''; q)n~k 

/gfc-n, 51-^-7(61626304), gl-"/(0l04),<?^-V(02a4),9^-V(«3a4) . . , 

1 ^1-V(?^ia4), gi-"/(&2a4), q^~^l{b,a,), q^-^^ / {aia^a-^a^) ' ' ' ^^•''^ 



Now substitute (j4.4p in (j4.7p and ()4.8p and let a; and y be given by ()4.5p . By similar computations 
as for obtaining (|4.6p we get as the limit of (14. 7p for g — )• 1 the following identity: 

^-^ Y" _ V" « + /^ + ^fe + ^ ^' p(-a-n-l,-/3-n-l), . p(",/3)/ n / . 

2 ) ~ ^^a + P + k + l + P + k + ^^^^^ ^"^^^ ^^'^^ 

Formula (14. 9p was earlier given by J. Koekoek & R. Koekoek [7, (21)]. As an alternative to 
their direct derivation (independently of the Askey- Wilson connection coefficients) one might 
also compute that 

!\{x - vY P^"''^^ (x) (1 - x)°(l + xf dx 

I-i{Pt''^\x))^ (1 - + x)l^ dx 

^ g + /3 + 2fc + 1 2" n! p(_„_„_i,_;3-n-i) 

a + l3 + k + l {a + P + k + 2) 



by substituting Rodrigues' formula for Jacobi polynomials in the numerator on the left-hand 
side, then performing repeated integration by parts, then using Euler's integral representation 
for hypergeometric functions and finally reversing the order of summation in the resulting ter- 
minating hyoergeometric series. 

From (j4.6p and ()4.9p we see (as also observed in [7J) that AB = I = BA, where A and B 
are the lower triangular matrices given for m > n > by 

_ {a + /3 + m+l)n p{Q+n,/3+n). ^ 

_ a + l3 + 2n + l m! <-a-m-i-i3-m-i), . (^-^O) 

a + p + n + l (a + /3 + n + 2)^ "^-^ ^^^^ 

In particular, we obtain from AB = I the identities (j2.15p and (j2.13p . while conversely from 
(I2.15P with (a,/3) running through all (a + j, /3 + j) (j £ Z>o) the full set of scalar identities in 
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AB = I for (a, /3) can be derived. Similarly we obtain from BA = I that 

g + /? + 2fc + 1 (g + ^ + l)fc p(«./3. N (-a-n-l,~/3-n-l) . . _ . s 

a + + l (a + /3 + n + 2)fc^'= '"-'^ (y) - Vo- (4.iij 

Formula (14. lip also follows from (14. 9|) by putting x = y, as already observed in [71 (22)]. 
Conversely (see [3 p. 13]), from (j4.1ip with {a, 13) running through all (a + j, /3 + j) (j G Z>o) 
the full set of scalar identities in BA = I for (a, /3) can be derived. 
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